Fundamental Parameters of He- Weak and He- Strong Stars * 



L. S. Cidale l , M. L. Arias 1 , A. F.Torres 1 , J. Zorec 2 , Y. Fremat 3 , A. Cruzado 1 



: 

o ■ 
o ■ 

(N ; 
a ■ 



^ . 

Or 
i < 

o ■ 

to: 



^1- 

q 

in 
o 
i> 
o 



1 Facultad de Ciencias Astronomicas y Geoffsicas, Uni- 
versidad Nacional de La Plata and Institute de Astroffsica 
de La Plata (CONICET), Paseo del Bosque S/N, 1900 La 
Plata, Buenos Aires, Argentina 

2 Institut d'Astrophysique de Paris, UMR7095 CNRS, 
Universite Pierre & Marie Curie, 98bis bd. Arago, 75014 
Paris, France 

3 Royal Observatory of Belgium, 3 Av. Circulaire, 1180 
Bruxelles, Belgium 

Abstract 

He-weak and He-strong stars are chemically peculiar AB 
objects whose He lines are anomalously weak or strong for 
their MK spectral type. The determination of fundamen- 
tal parameters for these stars is often more complex than 
for normal stars due to their abundance anomalies. 

We discuss the determination of fundamental parame- 
ters: effective temperature, surface gravity, and visual and 
bolometric absolute magnitudes of He- weak and He-strong 
stars. We compare our values with those derived inde- 
pendently from methods based on photometry and model 
fitting. 

We carried out low resolution spectroscopic observations 
in the wavelength range 3400-4700 A of 20 He-weak and 
8 He-strong stars to determine their fundamental parame- 
ters by means of the Divan-Chalonge-Barbier (BCD) spec- 
trophotometric system. This system is based on the mea- 
surement of the continuum energy distribution around the 
Balmer discontinuity (BD). For a few He- weak stars we 
also estimate the effective temperatures and the angular 
diameters by integrating absolute fluxes observed over a 
wide spectral range. Non-LTE model calculations are car- 
ried out to study the influence of the He/H abundance 
ratio on the emergent radiation of He-strong stars and on 
their T c q determination. 

We find that the effective temperatures, surface gravi- 
ties and bolometric absolute magnitudes of He- weak stars 
estimated with the BCD system and the integrated flux 
method are in good agreement between each other, and 
they also agree with previous determinations based on sev- 
eral different methods. The mean discrepancy between the 
visual absolute magnitudes derived using the Hipparcos 
parallaxes and the BCD values is on average ± 0.3 mag for 
He-weak stars, while it is ± 0.5 mag for He-strong stars. 



'Observations taken at CASLEO, operating under agreement of 
CONICET and the Universities of La Plata, Cordoba and San Juan, 
Argentina 



For He-strong stars, we note that the BCD calibration, 
based on stars in the solar environment, leads to overesti- 
mated values of T c g. By means of model atmosphere cal- 
culations with enhanced He/H abundance ratios we show 
that larger He/H ratios produce smaller BD which natu- 
rally explains the T e ff overestimation. We take advantage 
of these calculations to introduce a method to estimate 
the He/H abundance ratio in He-strong stars. The BD 
of HD 37479 suggests that the T e g of this star remains 
fairly constant as the star spectrum undergoes changes in 
the intensity of H and He absorption lines. Data for the 
He-strong star HD 66765 are reported for the first time. 

1 Introduction 

Ap/Bp stars are upper main sequence objects that show 
an abnormal abundance enhancement of some chemical 
species in their atmospheres as compared to "normal" A 
and B dwarfs of the same effective temperature (Jaschek 
& Jaschek [26]). Besides these chemically peculiar (CP) 
stars, which encompass classical A magnetic stars (CP1), 
Si, Cr and SrCrEu Ap stars (CP2) and HgMn stars (CP3), 
there are also CP B-type objects having abnormally weak 
or abnormally strong lines of Hci. They are called He- 
weak (CP4) and He-strong stars, respectively. The He- 
strong stars are considered the high-temperature extension 
of classical Ap/Bp stars (Osmer & Peterson |45j ) . 

Most CP stars have strong and variable dipolar or 
quadrupolar magnetic fields, whose interaction with the 
gravitational and radiative diffusion processes is assumed 
to induce non-homogeneous distributions of different 
chemical elements. However, in some cases, this mecha- 
nism alone cannot account for all abundance and isotopic 
anomalies (Profntt et al. [46]). 

Accurate values of fundamental parameters (effective 
temperature, surface gravity, visual and bolometric ab- 
solute magnitudes) are of crucial importance to study the 
evolutionary status of stars and the physical processes that 
take place in their atmospheres/interiors. They also en- 
able us to explore the appearance of their anomalies in 
connection with the physical characteristics of their local 
galactic environment (kinematics, magnetic fields, average 
content of chemical elements, etc.). In particular, these 
parameters are useful to study the evolution of magnetic 
fields in CP stars. 

Since the photosphcric abundance anomalies modify the 
flux distribution, the colours of CP stars can be affected 
by the concomitant blanketing effect. The use of standard 
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methods based on calibrations of colour indices in terms 
of T ff can then lead to erroneous results. Several authors 
have already discussed the difficulties of determining the 
effective temperature of these stars (e.g. Leckrone et al. 
[36] : Hauck & North [21]; Napiwotzki et al. |;43J ; Sokolov 
49J). From these studies it is apparent that the use of 
the Geneva or Stromgren photometry to determine T e g 
of CP stars is much more complex than for normal stars. 
Napiwotski et al. ([13]) have discussed different uvby/3 
photometric calibrations to derive T e g and logg. They 
recommended the use of Moon & Dworetsky's (1985) cal- 
ibration corrected for gravity deviations. 

Other alternative methods have also been proposed. For 
instance, the infrared flux method (IRFM) which mini- 
mizes the blanketing effect (Blackwell & Shallis [3]), mea- 
sures the ratio of the total integrated flux to a monochro- 
matic IR flux. This method is often used to determine 
angular diameters and effective temperatures of normal as 
well as of peculiar stars. Very frequently, the fundamental 
parameters are also determined by fitting synthetic spec- 
tra to the observed absolute energy distributions, or to line 
profiles, particularly the H/3 line. For CP2 stars, which are 
generally too blue for their spectral type, a supplementary 
correction to the effective temperature is needed (Stepieh 
& Dominiczak [ST]). 

Sokolov ([H5]) proposed the Balmer continuum slope 
near the Balmer jump as a tool to determine the T e g of 
CP stars. The method takes advantage of the fact that 
this slope is quite similar for normal and chemically pe- 
culiar main sequence stars. Although this method is easy 
to apply, it requires the stellar spectrophotometry to be 
corrected both for ISM reddening and for possible effects 
due to the abnormal energy distribution of CP stars. 

The surface gravity parameter log g is often determined 
through theoretical line profile fittings or by using evolu- 
tionary tracks. To enter the evolutionary tracks the ab- 
solute bolometric magnitude Mboi and T G ff of stars are 
needed. Most frequently, Mboi is determined with pho- 
tometric data and trigonometric parallaxes which provide 
the visual absolute magnitude that, in turn, is transformed 
into Mboi with the help of a bolometric correction calibra- 
tion (North et al. [H]; Gomez et al. [2D], Hubrig et al.[2"3). 
Kochukhov & Bagnulo [30]). 

Not only are these methods difficult to apply to CP stars 
but also none of them give simultaneously T Q g, logg, My 
and Mboi- Particularly, a) the photometry-based meth- 
ods need to take into account line blanketing effects, or 
they rely on corrections from chemical anomalies that are 
applied to colours; b) the measurement of the Balmer 
continuum slope is based on data previously corrected 
for energy redistributions and/or reddening effects; and 
c) the determination of visual absolute magnitudes from 
where the bolometric magnitudes are obtained, needs well- 
determined parallaxes and effective temperatures. 

Since most of these methods were applied to CP2 and 
CP3 categories, in the present contribution we focus our 
attention on He-weak and He-strong type stars. We note 



that information on He-strong stars is particularly scarce. 
In this work, we propose the use of the BCD spectrophoto- 
metric system (Barbier & Chalonge 2; ; Chalonge & Divan 
[7]) as an alternative tool to infer spectral types, effective 
temperatures, surface gravities, visual and bolometric ab- 
solute magnitudes of CP stars. This method has the ad- 
vantage of being based on direct measurable quantities of 
the continuum energy distribution in wavelengths near the 
the Balmer discontinuity (BD), which do not need cor- 
rections for interstellar and/or circumstellar extinctions 
and cannot be altered for beforehand unknown chemical 
anomalies. These quantities are strongly sensitive to the 
ionization balance of stellar atmospheres, so that they are 
excellent indicators of T e g and logg. Since these param- 
eters certainly depend somehow on the chemical proper- 
ties of stars, it seems useful to inquire whether or not we 
can also draw from them quantitative indications on their 
chemical anomalies. 

In §2 we present the observational material. In §3 we 
give a summary about the BCD spectrophotometric sys- 
tem and complementary techniques used to test against 
the BCD system. §4 is devoted to the determination 
of fundamental parameters of 20 He-weak stars with the 
BCD system. We also include our determinations of effec- 
tive temperatures, angular diameters and bolometric mag- 
nitudes done with the integrated flux method. We discuss 
the results with those obtained elsewhere independently. 
§5 gives fundamental parameters derived for 8 He-strong 
stars. Non-LTE model low resolution stellar spectra are 
also used to obtain the effective temperature and to ana- 
lyze the effect of enhanced He/H abundance ratio on the 
BD. We discuss the difficulties encountered in determin- 
ing the fundamental parameters of He-strong stars and the 
variation of HD 37479 spectrum mainly in view of its T e g 
estimate. §6 shows the evolutionary status of the stud- 
ied He-peculiar stars. Due to the variability observed in 
the He-strong group we add an appendix containing line 
equivalent widths. Conclusions are given in §7. 



2 Observations 

We obtained low resolution spectra of 28 CP stars: 20 
He- weak and 8 He-strong stars (see Table [TJ . Most stars 
were observed only once. Observations were carried out 
at the Complejo Astronomico El Leoncito (CASLEO), San 
Juan, Argentina, with the 2,15 m telescope and the Boiler 
& Chivens spectrograph on August 30-31, September 1-3 
and September 28-30, 2004, and on January 18-19, 2006. 
We used a 600 1mm- 1 grating (# 80), a slit width of 250 
/i and a CCD detector of 512 x 512 pixels. The obtained 
spectra cover the range 3400-4700 A with an effective res- 
olution of 4.53 A each 2 pixels. 

Bias, flat field, He-Ne-Ar comparison and spectropho- 
tometric flux standard star spectra were also secured to 
perform wavelength and flux calibrations. The reduction 
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of observations was made with the iraf0 software package 
(version 2.11.3) and all the spectra were corrected from at- 
mospheric extinction. 



3 The Methods 

3.1 The BCD spectre-photometric system 

The fundamental parameters of the program stars were 
mainly derived by using the BCD spectrophotometric sys- 
tem. This system is based on the observable parameters 
(Ai, D). D measures the Balmer jump at A3700 A and it 
is a strong indicator of the effective temperature. Ai gives 
the mean spectral position of the Balmer jump and it is 
related to the surface gravity. 

The BCD spectrophotometric system has several ad- 
vantages: 1) it is not affected either by interstellar or cir- 
cumstellar extinction (Zorec & Briot [BDJ), 2) it is based 
on measurable parameters of the visual continuum energy 
distribution near the BD which are thus related, on aver- 
age, to physical properties of deeper photospheric layers 
than any classification system based on spectral line mea- 
surements, 3) it is possible to derive information on most 
of the fundamental parameters, 4) this method is easy and 
direct to apply, and 5) it can be applied with high accu- 
racy to normal and many peculiar classes of B-type stars. 
Zorec & Briot ([BD]), Cidale et al. Q10J) and Zorec et al. 
([55]) applied it successfully to Be and B[e] objects, even 
though these objects may display a second BD that can 
be either in emission or in absorption (Divan [13 ). 

To determine the parameter D — log(F +3 7 o/-F-37oo) we 
normalize the observed energy distribution with a Planck- 
ian function, which rectify the energy distribution at both 
sides of the BD, and then, we extrapolate the Paschen con- 
tinuum to obtain the flux at ^+3700- To determine i ? _37oo, 
we use the flux level where the higher members of Balmer 
lines merge together. The use of Balmer lines to determine 
the low level of the Balmer jump has the advantage of giv- 
ing an estimate of the BD which is not affected by possible 
variations of the intensity of the Balmer continuum. 

We have measured D and Ai for the program stars and 
determined their spectral type, T c g, logg, and absolute 
magnitudes My and Mboi (see details in §4 and §5) using 
the respective calibrations of (Ai,D) by Divan & Zorec 
(El), Zorec (P|), and Zorec & Briot ([BD]). The uncer- 
tainties on T e ff, loggf, My and -Mboi are mainly related to 
the errors of the BCD observational quantities (Xi,D). 
Typical values are a (D) ~ 0.015 dex and <r(Ai) ~ 5 A, 
which produce AT e g ~ ±500 K for the late B-type stars 
and ±1500 K for the early B-types (T eff > 20000 K), 
A log g ~ ±0.2 dex, and AM ~ ±0.3 mag both for My 
and Mboi- 

1 IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for Re- 
search in Astronomy (AURA), Inc., under cooperative agreement 
with the National Science Foundation 



Table 1: Program He- weak (top) and He-strong (botton) 
stars. In the first two columns we give the star identifiers. 
Column 3 shows the date of observation. The peculiarity 
class, and the average quadratic effective magnetic field 
are presented in columns 4 and 5 and were taken from 
Bychkov et al. (2003) 



HD 


HR 


Date 


Peculiarity 


<B e > 






dd/mm/yy 




Gauss 


5737 


280 


01/09/04 


SrTi He-w 


324.0 






03/09/04 










04/09/04 






19400 


939 


01/09/04 


He-w 


206.7 


22470 


1100 


01/09/04 


Si He-w 


732.9 


22920 


1121 


30/09/04 


Si4200 He-w 


307.1 






01/10/04 






23408 


1149 


03/09/04 


Hg He-w 


410.0 


28843 


1441 


03/09/04 


Si He-w 


344.2 


49333 


2509 


03/09/04 


Si He-w 


618.4 






01/10/04 






49606 


2519 


01/10/04 


MnHgSi He-w 


916.0 


51688 


2605 


01/10/04 


Hg He-w 


550.2 


57219 


2790 


04/09/04 




- 






18/01/06 






74196 


3448 


18/01/06 


He-w 


- 


125823 


5378 


18/01/06 


Si He-w 


469.3 


142301 


5912 


01/10/04 


Si He-w 


2103.6 


142990 


5942 


02/09/04 


He-w 


1304.3 






30/09/04 






144334 


5988 


01/10/04 


Si He-w 


783.2 


144661 


5998 


01/10/04 


He-w 


542.0 


144844 


6003 


01/10/04 


Si He-w 


318.1 


162374 


6647 


01/10/04 


He-w 


269.8 


175362 


7129 


03/09/04 


SiMn He-w 


3569.9 


202671 


8137 


01/09/04 


SrTi He-wm 


183.0 


37017 


1890 


01/09/04 


He-s 


1488.1 






30/09/04 






37479 


1932 


02/09/04 


He-s 


1907.9 






04/09/04 






37776 


BD-01 1005 


04/09/04 


He-s 


1259.7 


58260 


CD-36 3578 


18/01/06 


He-s 


2291.2 


60344 


CD-23 5673 


18/01/06 


He-s 


334.9 


64740 


3089 


03/09/04 


He-s 


571.9 






30/09/04 






66765 


CD-47 3537 


18/01/06 


He-s 




264111 


BD+04 1447 


18/01/06 


He-s 
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3.2 The integrated-flux method 

In addition to the BCD classification system we used the 
integrated-flux method to derive fundamental parameters 
for those He-weak stars of our program with absolute 
fluxes observed over a wide spectral range. We perform 
a determination of the effective temperature, angular di- 
ameter (0), and bolometric magnitude by means of the 
relation: 

F hoi = {6 2 /A)a^T^ (1) 

where Fboi is the absolute bolometric flux reduced to the 
distance to the Earth and corrected for interstellar extinc- 
tion; crg B is the Stefan-Boltzmann constant. Relation ([T]) 
is iterated starting with an approximate T G ff until the dif- 
ference between two succesive iterations becomes AT e g ~ 1 
K. This implies from 10 to 30 iteration steps. At each 
step, 8 is calculated using the F\/F mo d e \ ratio of fluxes in 
the wavelength interval AA7000 - 8000 A. The F x /F nwdcl 
ratios are fairly insensitive to log<?. The observed fluxes 
F\ are corrected for ISM extinction and the i^ mo dei corre- 
sponds to the running T e ff . At each step we also determine 
the amount of non-observed far-UV and far-IR energy as 
a function of the T e g under iteration using LTE models 
(Kurucz J33J), which we add to the integrated observed 
flux to obtain Fboi- The wavelength interval of observed 
absolute fluxes ranges from 1300 A to 1 fim. The fluxes 
are observed in the far-UV by the TD1 and IUE satel- 
lites, and from visible to near-IR with the 13-colour narrow 
band photometry calibrated in absolute fluxes (Johnson & 
Mitchell [1H]). 

3.3 Model fitting 

In the case of He-strong stars we also determine T e ff, 
and log<? by computing NLTE synthetic spectra with 
the TLUSTY and SYNSPEC computer codes (Hubeny 
& Lanz [33], and the references therein) assuming H-He 
model atmospheres and He/H ratios of 0.1, 0.2, 0.5 and 
1.0 . The atomic models we used are basically those pro- 
vided on the TLUSTY website for Hi (9 levels), Hei (20 
individual levels) and Hen (20 levels). The Stark broad- 
ening of the He I transitions close to the (ls2p) 3Po ion- 
ization limit is taken from Dimitrijevic & Sahal-Brcchot 
([llj.[T2j). or it is computed using an approximate relation 
proposed by Freudenstein & Cooper ([16]). In all cases, 
a microturbulent velocity of 2 km s _1 is assumed. To 
have synthetic spectra comparable with the ones observed 
at CASLEO, we resampled them to the resolution of the 
observed energy distributions. 

4 He-weak stars 

Table [2] gives measurements and fundamental parameters 
derived from the BCD spectrophotometric system. Ta- 
ble [3] reproduces the values of the fundamental parameters 
derived by using integrated absolute fluxes. 



Table 2: He- weak stars: measured (Ai,£>) BCD quanti- 
ties and the corresponding fundamental parameters: T e ff, 
\ogg, spectral type, My, and Mboi determined using the 
calibrations of (Ai, D) for OB stars of solar He/H ratio. 



HD 


Ai 


D 


Toff 


log 


Sp.T 


My 


Mboi 




A 


dex 


K 


dex 




mag 

o 




5737 


OO T A 


0.281 


1 OTA A 

12700 


O A C 

2.95 


B6III 


A C 

—2.5 


—3.4 


~i f\ a aa 

19400 


a o on 
48.30 


A OOA 
0.321) 


i oaaa 
13000 


O AA 

3.90 


B6-7V 


A O 

— 0.8 


1 C 

— 1.5 


OO A TA 

224/0 


58.68 


a one 
0.296 


i a aaa 
14000 


A 1 

4.21 


B6V 


A Q 

— 0.3 


— 1.6 


22920 


AC 1 A 

46.10 


A 077 

0.277 


1 A AAA 

14200 


O OA 

3.80 


B6IV 


— 1.5 


o c 

—2.5 


OO A AO 

23408 


on ao 

39.02 


A A AO 

0.408 


1 1 A C A 

11250 


O OA 

3.20 


B9III 


— 1.1 


— 1.6 


O O O A O 

28843 


50.67 


0.273 


14700 


3.95 


B5IV 


— 1.1 


—2.2 


49666 


c 1 oa 

61.20 


A OCA 

0.250 


1 C OAA 

15800 


A OO 

4.23 


B5V 


A T 

—0.7 


O O 

—2.3 


49606 


Ad CO 

46. oo 


A O A A 

0.340 


IOC A A 

13500 


O OA 

3.8U 


B6IV 


1 o 
— 1.2 


O A 

— 2.0 


ri poo 

51688 


AC\ AO 

49.9(5 


A O A A 

0.340 


IOC A A 

12500 


O AC 

3.95 


B7-8IV 


A £2 

—0.6 


1 o 

— 1.2 


c to i n 

57219 


C 1 A 1 

51.91 


A 1 AO 

0.193 


IOC A A 

18500 


O AO 

3.92 


B2.5IV 


O A 

—2.0 


O T 

—3.7 


74196 


54.78 


0.299 


13750 


4.13 


B6V 


—0.5 


— 1.6 


1 o r* o o o 

125823 


57.64 


0.184 


19400 


4.10 


B2V 


— 1.7 


—3.6 


142301 


65.60 


0.252 


15700 


4.30 


B5V 


-0.7 


-2.2 


142990 


67.21 


0.208 


18500 


4.27 


B3V 


-1.2 


-3.0 


144334 


63.76 


0.280 


14500 


4.30 


B5-6V 


-0.4 


-1.8 


144661 


60.31 


0.273 


15000 


4.23 


B5V 


-0.5 


-2.0 


144844 


65.60 


0.346 


12100 


4.32 


B8V 


+0.4 


-0.8 


162374 


51.11 


0.247 


15700 


3.95 


B5V 


-1.3 


-2.7 


175362 


64.75 


0.219 


17500 


4.25 


B3V 


-1.0 


-3.0 


202671 


44.37 


0.303 


13200 


3.70 


B6IV 


-1.3 


-2.0 



Note: The Ai parameter is given in A — 3700 A 



HD 125823 is the only star whose Hei lines vary 
in strength between anomalously weak and anomalously 
strong. It was included in Table [2] because the star pre- 
sented a He- weak phase at the date of our observations. 

4.1 T c g determination 

In order to test the accuracy of the BCD effective tem- 
peratures a comparison with the values obtained by other 
methods is performed. Table|4]shows the effective temper- 
atures with their corresponding uncertainties (when avail- 
able) derived from the BCD system and from the litera- 
ture. 

Except for few stars, the effective temperatures derived 
with the BCD method have no significant differences with 
the values derived from Geneva and other multicolour pho- 
tometric systems, the IRFM and LTE line-blanketed mod- 
els. The mean discrepancies found between our data and 
those derived by other authors are lower than 800 K. The 
best agreement is obtained with the IRFM (Lanz [34]). 
the mean temperature difference being |A T e g \ ~ 500 K. 
We have also obtained a good agreement between our own 
determinations of T c s{IF) by means of the integrated-flux 
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Table 3: He-weak stars: our determinations of fundamen- 
tal parameters: T c s, 6,Ft, i, -^boi and the E(B — V) color 
excess used for the ISM extinction correction applying the 
integrated-flux method. 



HD Teff(IF) #(rad) F bol M bol E B _ V 

K xlO~ 10 xlO -7 mag mag 

5737 12490 ±143 17.6 ±0.2 10.6 ±0.6 -3.1±0.4 0.00 

57219 18286 ± 125 8.9 ± 0.1 12.5 ± 0.5 -3.8±0.3 0.05 

125823 17594 ± 172 12.8 ± 0.2 22.2 ± 1.7 -2.9±0.2 0.01 

142301 15084 ± 84 8.8 ±0.1 5.7 ± 0.2 -1.6±0.4 0.12 

144334 15099 ± 77 8.1 ±0.1 4.8 ± 0.2 -1.6±0.3 0.10 

144661 13421 ± 121 7.6 ± 0.5 2.7 ± 0.4 -0.4±0.2 0.09 

144844 12448 ± 50 10.1 ±0.1 3.5 ± 0.1 -0.9±0.2 0.12 



Table 4: Comparison of the BCD effective temperatures 
with those determined by other authors for the program 
He-weak stars. Column 1 gives the HD number; column 
2, our T c fr(BCD) values; column 3, T e g(G) (Glagolevskij 
[18J) determined from the parameters Q in the UBV sys- 
tem and X in the multicolour system; column 4, T ff (KB) 
determined by Kochukhov & Bagnulo ([10]) using the cal- 
ibration of the Geneva photometry (Golay [H5]) and the 
method suggested by Hauck & North ([21]); column 5: 
T cff (HG) from Hunger & Groote ([25]) using the IRFM, 
column 6: T e ff (L) obtained by Lanz ([M]) using the IRFM, 
column 7: T cff (HN) from Hauck & North ([ST]) using 
the Geneva photometry, and column 8, other values of 
T e ff ± AT found in the literature, most of them deter- 
mined by model fitting 



HD 



T ctf x 10 3 K 



BCD 



KB 



HG 



method and the BCD system. Comparing Tables [SJ and [3] 
we found a mean value of AT e g =±700i"T. 

Figure [TJ displays the comparison of our T c ff(BCD) val- 
ues vs. T e ff(G) (crosses), T e ff(KB) (open circles), and 
T c ff(IF) (open triangles). In the same figure, it is also 
plotted a reference line with slope 1. We selected the 
mentioned works to perform the comparison, since the 
number of stars in common is large. The good correla- 
tion found between the BCD system and other methods 
confirms that the BD of He-weak stars behaves like that 
observed in normal B stars. For these stars we can then 
safely use the existent calibrations, done for He/H w 0.1, 
to determine their fundamental parameters. 

From our results we conclude that the behaviour of the 
Balmer jump in He-weak stars is different from that of Ap 
stars (CP2 group). The photometric anomalies of Ap stars 
(Gerbaldi et al. [T7]; Hauck & North [SSJ) may introduce 
noticeable differences between the BCD and MK spectral 
classifications (Chalonge & Divan [9]). 



4.2 The surface gravity 

Surface gravity determinations for He-peculiar stars are 
hardly common. Hunger & Groote (1999) determined 
log g values from hipparcos parallaxes together with 
tracks from stellar evolution, Leone & Manfre (1996) es- 
timated values using synthetic H line profiles while Leone 
et al. (1997) used Moon & Dworetsky's (1985) calibration 
and applied the corrections suggested by Napiwotzki et al 
(1993). Table [5j compares, for each object, the log g val- 
ues determined by means of the BCD system with those 
derived by the above mentioned works. We find an excel- 
lent agreement with the values given by Hunger & Groote 
(1999), | A log g\ < 0.13 dex, for 11 stars in common and a 
good agreement with Leone's works, |Alog g\ = 0.23 dex 
for six stars in common. 



5737 12.7±0.5 13.50±0.3 13.2±0.4 14.4±0.1 



19400 13.0±0.5 12.90±0.3 



22470 
22920 
23408 
28843 
49333 
49606 
51688 
57219 
74196 
125823 
142301 
142990 
144334 
144661 
144844 
162374 
175362 
202671 



14.0±0.5 
14.2±0.5 
11.2±0.5 
14.7±0.5 
15.8±0.5 
13.5±0.5 
12.5±0.5 
18.5±0.5 
13.7±0.5 
19.4±0.5 
15.7±0.5 
18.5±0.5 
14.5±0.5 
15.0±0.5 
12.1±0.5 
15.7±0.5 
17.5±0.5 
13.2±0.5 



13.45±0.3 
14.40±0.3 
12.30±0.3 
14.53±0.3 
16.60±0.3 
13.80±0.3 
13.05±0.3 
16.35±0.3 
13.30±0.3 
19.53±0.3 
16.50±0.3 
17.80±0.3 
15.40±0.3 
15.00±0.3 
12.30±0.3 
17.30±0.3 
17.00±0.3 
14.10±0.3 



13.0±0.4 
13.9±0.4 
12.0±0.4 
13.9±0.4 
16.4±0.5 



17.7±0.5 
15.6±0.4 
16.5±0.5 
14.7±0.4 



13.5±0.4 
14.3±0.3 



15.0±0.5 
15.7±0.4 



14.2±0.6 
18.4±0.4 
15.9±0.3 
16.9±0.5 
15.2±0.5 



HX 



13.79±0.14 14.25 



14.80±0.19 
15.94±0.22 



13.80 
15.28 
16.81 



Others 

13.6±0.2° 
14.0 6 

13.3±0.5" 
13. 9 C 

13.4±0.5° 

14.1±0.5 a 



14.8° 



17.1° 



19.60 
17.26 



19.01±0.23 
15.99±0.15 
18.02±0.23 
14.58±0.23 16.25 



16.1±0.3 17.84±0.22 17.37 
16.5±0.4 16.38±0.20 18.24 



16. 5 e 
13.1±0.2" 



a: Lconc & Manfre ([37]) or Leone et al.([38j), b: Cayrel De Strobel et al.([6]) 
c: Jaschek & Jaschek (1987b), d: Kroll (SFTj), e: Mcgessicr (1988b) 



4.3 Visual and bolometric absolute mag- 
nitudes 

Average absolute magnitudes for Ap/Bp stars were ob- 
tained by North et al. (1997) and Megessier (1988a), 
and for individual objects were reported by Gomez et al. 
(1998) and Kochukhov & Bagnulo (2006). Using photo- 
metric measurements and HIPPARCOS astrometric data, 
Gomez et al. (1998) obtained that the He-weak stellar 
group has My —0.2 mag and that the intrinsic dis- 
persion varies from 0.6 to 0.8 mag. All these stars are 
clearly distributed in the Gould Belt and their kinematic 
behaviour is similar to that of the non-peculiar thin disk 
stars younger than about 10 9 years. 

In our sample of 20 He- weak stars, 16 of them have abso- 
lute magnitudes determined by the LM algorithm (Gomez 
et al. [20]), 10 are in common with Glagolevskij ([18]) and 
11 are in common with Kochukhov & Bagnulo's ( [30] ) 
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Figure 1: Comparison of the effective temperatures de- 
rived from the BCD system with those obtained photo- 
metrically by Glagolevskij ([18]), Kochukhov & Bagnulo 
( 30J) and in this work with the integrated flux method 
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Figure 2: Comparison of the visual absolute magnitudes 
derived from the BCD system with those derived by 
Gomez et al. (1998), Glagolevskij (2002) and Kochukhov 
& Bagnulo (2006). He- weak stars are represented by open 
symbols and He-strong stars by filled symbols 



Table 5: He-weak stars: comparison of BCD log g values 
with those obtained by other authors 



HD 


logg(BCD) 


logg(HG) 


logfif(Le) 




dex 


dex 


dex 


5737 


2.95±0.2 


3-48t°;ii 


3.20±0.10 


19400 


3.90±0.2 




3.76±0.25 


22470 


4.21±0.2 


o Q7+0.09 
o i " -0.08 


4.15±0.25 


22920 


3.80±0.2 


O O + 0.13 
— U. 1 1 


3.72±0.25 


23408 


3.20±0.2 







28843 


3.95±0.2 


4.33+Hs 


4.25±0.25 


49333 


4.23±0.2 


4.18+"-^ 

— U. 1U 


4.08±0.25 


49606 


3.80±0.2 


3.89±0.25 


51688 


3.95±0.2 


— 


— 


57219 


3.92±0.2 


— 


— 


74196 


4.13±0.2 


4.14lg"» 


— 


125823 


4.10±0.2 


4.16± ) : ° 1 8 7 


_ 


142301 


4.30±0.2 


4.29 + ° : i? 


_ 


142990 


4.27±0.2 


4.20 +0 j : ° 1 


_ 


144334 


4.30±0.2 


4.29+° : ° 


_ 


144661 


4.23±0.2 




- 


144844 


4.32±0.2 






162374 


3.95±0.2 


°- JO -0.14 




175362 


4.25±0.2 




3.70±0.10 


202671 


3.70±0.2 




3.40±0.10 



Column 2: this work, Column 3: Hunger & 
Groote (1999), Column 4: Leone & Manfre ([57]) 
or Leone et al.([55]) 



work (see Table E]) ■ Our determinations of visual abso- 
lute magnitudes are in a very good agreement with those 
obtained by the authors mentioned above (see Figure [2]). 
The mean discrepancy between the BCD method and the 
LM algorithm for the He-weak stars is 0.26 mag. We also 
found a good agreement with My obtained by Kochukhov 
& Bagnulo (2006), AM V = ± 0.3. Instead, the mean 
discrepancy between the bolomctric magnitudes derived 
from the BCD system and the values given by Glagolevskij 
([18J), using the parameter (3 of the multicolour photome- 
try, is slightly larger ~ 0.44 mag. This discrepancy may be 
related to the bolometric correction, since when transform- 
ing the bolometric magnitudes obtained by Glagolevskij 
(|18j) into visual absolute magnitudes (which are given in 
Table [5]), the discrepancy become smaller. 

Comparing the bolometric absolute magnitudes ob- 
tained by us with the BCD system and the integrated- 
flux method, we find a mean discrepancy AAiboi = ± 0.5 
mag. However, a noticeable disagreement appears for HD 
144661 and for the He-variable HD 125823. The large dif- 
ference in the Mboi observed for HD 144661 cannot be 
recovered by simply changing the color excess E(B — V); 
by doubling the E(B — V) the bolometric luminosity would 
increase AL/ L ^ 0.5, while we need AL/L^ZA to make 
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Table 6: He-weak stars: Comparison of the BCD absolute 
visual magnitudes with those derived by other authors. 
The uncertainties of the values are quoted when they are 
available 



HD M v (BCD) M V (G) M V (KB) M v (Go) 

2.0 
0.7 
0.0 
1.3 
1.3 
0.0 
0.5 
1.0 
0.8 

0.4 
1.4 
0.5 

0.5 



1.4 
0.5 
1.1 

Column 2: this work, Column 3: Glagolevskij (2002) 
Column 4: Kochukhov & Bagnulo (2006), and 
Column 5: Gomez ct al. (1998) 



both estimates agree. Therefore, we may suggest that 
some undetected flux fading is taking place. This star 
deserves perhaps further thorough analysis. 

5 He-strong stars 

The shape of the Balmer continuum of He-strong stars 
is somewhat different from that predicted by models with 
standard solar helium abundances. The He-strong spectra 
exhibit intense He I lines and a remarkable convergence of 
the last members of the series to the bound free discon- 
tinuity at 3422 A, which corresponds to the ionization of 
the 2-3 Po level (see Figs. [3] and d]). The He-strong stars 
deserve a special discussion. They may display Ha emis- 
sion and spectroscopic and photometric variability. Not 
only do they show variations in line intensities, radial ve- 
locities, luminosity and colours but also in the magnetic 
field strength (Renson et al. [4T]). 

We have attempted at determining fundamental param- 
eters of these objects by means of the BCD system. How- 
ever, due to the strong He I lines, a precise determination 
of the Paschen continuum level becomes difficult. Never- 
theless, the continuum fluxes at AA4055 and 4566 A are un- 



Table 7: He-strong stars: measurements and fundamen- 
tal parameters derived from the BCD spectrophotometric 
system 



HD 


Ai 


D 




logff 


Sp.T 


My 


M bol 




A 


dex 


K 


dex 




mag 


mag 


37017 


68.72 


0.145 


21700 


4.22 


B1.5V 


-2.0 


-4.5 


37479 


87.77 


0.107 


25300 


4.30 


B0.5V 


-2.0 


-5.3 


37776 


74.68 


0.107 


25400 


4.25 


B0.5V 


-2.2 


-5.5 


58260 


66.05 


0.149 


21400 


4.20 


B2V 


-1.9 


-4.4 


60344 


52.97 


0.111 


24600 


3.80 


B1IV 


-3.4 


-5.7 


64740 


77.97 


0.115 


24500 


4.23 


B1V 


-2.1 


-5.2 


66765 


60.79 


0.160 


20200 


4.11 


B2V 


-2.0 


-4.4 


264111 


54.78 


0.092 


26700 


3.85 


BOIV 


-3.7 


-6.5 



Note: The Ai parameters are given in A — 3700 A 



Table 8: Stellar parameters for the reference He-strong 
stars obtained with NLTE model atmospheres. We list 
the HD numbers (column 1), the adopted T c g, log<?, and 
Hc/H ratio (columns 2, 3, and 4, respectively), and V^sini 
(column 5). Vsmi values are taken from Veto et al. (|53j). 

Wal born ([55]), and Zboril et al. ([57]) 

"HD T cff (NLTE) log^ Uc/U Vsmi 





[K] 


[dex] 




km s 1 


37017 


21000 


4.50 


0.5 


95 


37479 


23000 


4.00 


0.5 


162 


37776 


23000 


4.00 


0.5 


95 


58260 


19000 


4.00 


1.0 


45 


60344 


24000 


4.50 


1.0 


55 


64740 


23000 


4.00 


0.5 


160 


66765 


20000 


4.00 


0.1 


100 


264111 


21000 


4.50 


1.0 


140 



affected by strong line absorptions (Adelman & Pyper pQ). 
Thus, we use these points, together with the last mem- 
bers of the Balmer line series as references to determine 
the Balmer jump. Table[7|gives our (Ai, D) measurements 
and the corresponding fundamental parameters. Data for 
the He-strong star HD 66765 is scarce and spectrophoto- 
metric measurements are published for the first time. 

The energy distributions of the observed He-strong stars 
were fitted with the TLUSTY and SYNSPEC computing 
codes to have another independent effective temperature 
estimate. All models were convolved with a rotational 
broadening function for the respective stellar Vsmi pa- 
rameters. In Figures [3] and [5J we show the best fit ob- 
tained for each studied star. These fits are for the (T e g, 
logg, H/He ratio) sets of parameters given in Table [S] 
The uncertainties on the parameters obtained by model 
fitting are on average: 1000 K for T e g and 0.5 dex in \ogg. 
However, for each He-strong star, the model that best fits 
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Figure 3: Energy distributions of He-strong stars, log .Fa, around the BD (dotted line), together with the best theo- 
retical fitting (solid line) . All spectra have the same flux scale. The convergence of the last He I line members towards 
3422 A clearly indicates the He I bound-free discontinuity 



the continuum energy distribution does not correctly fit 
the intensity of H lines. A better fitting of H lines would 
require models with larger atmospheric temperature gradi- 
ents than the ones used in NLTE blanketed plane-parallel, 
in hydrostatic and radiative equilibrium, classical models. 
Such differences in the atmospheric temperature structure 
could be due not only to possible non-uniform distribu- 
tions of the abnormally abundant chemical elements in 
latitude and/or in depth, but also to a non-negligible role 
of the magnetic field in cooling the upper atmospheric lay- 
ers. 

5.1 T c g determination 

We have sought for previous determinations of effective 
temperatures in the literature. Most of the T c ff values were 



in works that deal with only a reduced number of stars and 
use different techniques (Hunger & Groote [35]; Leone et 
al. [35]; Zboril et al. [55]; Glagolevskij [TSj; Theodossiou 
& Danezis [52] ; Adelman & Pyper [I] ; Kaufmann & Theil 
[29]). 

Table [H] compares the BCD effective temperatures with 
those obtained by Glagolevskij ([H]), Zboril et al. f[55]) 
and Kaufmann & Theil ([29 ). In this Table, the non- 
BCD temperature determinations are from theoretical 
LTE model-fitting (Kurucz [33]), or from the Geneva 
photometry. We have chosen only these works because 
they treat a non negligible number of stars that are in 
common with our program stars. We found large dis- 
crepancies among the T e g estimated by different authors: 
2000~AT Gff -5000 K. In the same way, our BCD deter- 
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Figure 4: Energy distributions of He-strong stars, log .Fa, around the BD (dotted line), together with the best theo- 
retical fitting (solid line) . All spectra have the same flux scale. The convergence of the last He I line members towards 
3422 A clearly indicates the He I bound-free discontinuity 



minations do not escape the rule. Moreover, we observed 
large differencies between the T e g derived with the BCD 
method and the NLTE models we used. In most cases, 
the T c ff(BCD) are systematically larger than the average 
Teg from other methods. This difference deserves a more 
detailed discussion that will be done using the stellar at- 
mosphere models in H5.4I 

5.2 The surface gravity 

As it is mentioned in H4.21 the log g values found in the lit- 
erature were derived by different methods. In Table [TU1 we 
show BCD log g values together with the determinations 
given by other authors. The resulting mean discrepancies 
in log g are 0.36, 0.21, 0.37 and 0.22 dex, respectively. The 
values of log g quoted in the atlas of Kaufmann & Theil 



(1980) are taken from different sources and their uncer- 
tainties are not accessible. 

The best agreement is obtained with values derived by 
Hunger & Groote (1999) from HIPPARCOS parallaxes to- 
gether with tracks from stellar evolution. 

5.3 Visual and bolometric absolute mag- 
nitudes 

The He-strong group typically has My ~ — 1.6 mag with a 
dispersion of about 1.2 mag (Gomez et al 1998) . 

In our sample of 8 He-strong stars, 6 of them have abso- 
lute magnitudes determined by the LM algorithm (Gomez 
et al. [20]), 5 are in common with Glagolevskij ([18]) and 
1 with Kochukhov & Bagnulo's ([30]) work (see Table fTTj) . 
In Table QTJ uncertainties for each value are listed when 
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Table 9: Comparison of the BCD effective temperatures 
with those determined by other authors for the program 
He-strong stars 



Table 11: Absolute visual magnitudes for the program He- 
strong stars 



HD 




T cp x 1 

± cfi A iu 


3 K 




BCD 


G 


Z 


K 


3701 7 

JIU1I 


21.7±1.5 


on i n_i_n q 

Z*U. 1U- 1— u . o 


19.2 


21.0 


^747Q 


OK o_|_1 C 

z<j . onz l . o 


99 4^,-1-0 3 

^ . 4: ZC U . O 


99 9 
_ _ . — 


9^! ^ 
zo . 


O i i i u 






91 8 

Z _L . O 


9fi n 

zu . u 


58260 


21.4±1.5 


19.70±0.3 


19.0 




60344 


24.6±1.5 


21.90±0.3 


21.7 


25. Oil. 


64740 


24.5±1.5 


23.85±0.3 


22.7 


23.5 


66765 


20.2±1.5 








264111 


26.7±1.5 


21.60±0.3 


23.2 


20.0 


Column 2: 


BCD syste 


m; Column 3: 





Glagolevskij (2002); Column 4: Zboril ct al. (1997); 
Column 5: Kaufmann &Theil (1980) 



Table 10: 
He-strong 

HP 

37017 
37479 
37776 
58260 
60344 
64740 
66765 
264111 



Comparison of surface gravity estimates for the 
stars 

logg(BCD) logg(Z) logg(HG) logg(LCM) logg(K) = 
A oo+n n A A^4-n in 4 12+ n34 a m+n or; A a 



4.22±0.2 
4.30±0.2 
4.25±0.2 
4.20±0.2 
3.80±0.2 
4.23±0.2 
4.11±0.2 
3.85±0.2 



4.45±0.10 
4.53±0.17 
4.52±0.07 
4.02±0.18 
4.48±0.18 
4.50±0.06 

4.54±0.06 



4.16 
3.55 



4.10 



-0.23 

+0.30 
0.30 
+0.47 
0.27 

+0.09 



4.02±0.25 
3.86±0.25 
4.11±0.25 
3.52±0.25 
3.50±0.25 
3.80±0.25 



4.4 
4.1 
4.0 

3.62 
3.9 

4.0 



Column 2: this work, Column 3: Zboril ct al. (1997), 

Column 4: Hunger & Groote (1999), Column 5: 

Leone et al. (1997), Column 6: Kaufmann & Thcil (1980) 



available. The values corresponding to My(G) were de- 
rived from Mboi given by Glagolevskij (2002) adding the 
respective bolometric correction. Our determinations of 
visual magnitudes are in good agreement with those ob- 
tained by the authors mentioned above (see Figure^ filled 
symbols), with the exception of HD 58260. We also found 
that for HD 64740 M V (BCD) differs significantly with the 
value obtained by Glagolevskij (2002). We find an aver- 
age discrepancy for all the stars in common of 0.53 mag 
between the My obtained with the BCD method and the 
LM algorithm. Comparing the values of Mboi(BCD) with 
the absolute bolometric magnitudes given by Glagolevskij 
([18j) we find an average discrepancy of 0.71 mag. This 
large difference may probably due to the bolometric cor- 
rection used by Glagolevskij ([H]) as we mentioned in M4.31 

5.4 Other relevant results 

a) - Synthetic energy distributions as a function of the 
He/H abundance ratio 

Since BCD determinations of T c g for the He-strong stars 
are systematically higher than those obtained by other 



HD 


Afy(BCD) 


M V (G) 


M v (KB) 


M v (Go) 


3701 7 


o n-i-n ^ 






9 9 
z . z 


3747Q 










3777^ 


9 9_i_n q 


9 7 
Z . t 




9 9 
z . z 


58960 


_i q_i_f) q 


-4.2 




-3 3 


60344 


-3.4±0.3 


-3.6 




-2.7 


64740 


-2.1±0.3 


-3.3 


-2.15±0.25 


-2.5 


66765 


-2.0±0.3 








264111 


-3.7±0.3 






-3.1 



Columns 2: this work, Column 3: Glagolevskij (2002), 
Column 4: Kochukhov & Bagnulo (2006), and 
Column 5: Gomez et al. (1998) 



authors, we calculate low-resolution synthetic spectra, in 
order to quantify the influence of the He/H abundance ra- 
tio on the emitted visual energy distribution. From these 
spectra we obtain the BCD parameters following the same 
procedure as for the empirical ones. Models are computed 
for effective temperatures close to those inferred empiri- 
cally for He-strong stars and only for the wavelength in- 
terval of the observed low-resolution spectra. 

An illustration of the effects produced on the visible 
energy distributions resulting from increased abundance 
ratios He/H is shown in Figure [5] The spectra are slightly 
displaced in order to see clearly the changes in the BD, 
where it is apparent its shortening as He/H is higher. In 
this figure we can also notice the changes on the visible 
Balmer energy distribution induced by enhanced He free- 
bound transitions and the presence of strong He I line ab- 
sorptions. In Table Q21 the obtained D values are given 
as a function of the model nominal (T c ff,logg, He/H) pa- 
rameters. Thus, increased He/H ratios produce smaller 
values of D. This reduction has only a second order de- 
pendence with log g. The lowering <5D of the BD, aver- 
aged over log g, as function of the effective temperature 
is given within 99.8% of approximation by the following 
interpolation relation: 



<5D = -0.056 



1 - 0.233 



To 1 



0.974 



(He/H) (dex) 



(2) 

Relation $2$ implies that the higher the T c g the smaller the 
<5D. However, as D oc 1.0/T e ff, small changes in D carry 
big differences in T c g. Since the BCD calibrations used to 
estimate the effective temperatures of He-strong stars are 
actually made for objects with He/H s=s 0.1, from Table [T2l 
and relation it becomes clear that the obtained BCD 
T e ff are overestimated. 

Models also predict a slight brightening of stars as the 
Hc/H ratio increases. Table [TBI reproduces the changes of 
the monochromatic magnitude at A4250 A with respect 
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Table 12: Model D parameters as a function of T e g, logg 
and He/H abundance ratio 



3800 4000 

Wavelength [A] 



Figure 5: Synthetic spectra displaying line and continuum 
changes according to increased He abundances. The spec- 
tra are slightly displaced. The NLTE model computations 
were performed for T e g = 23000 K, log g = 4.0 dex and 
ratios H/He= 0.1 (standard), 0.5 and 1.0 



to models for He/H = 0.1 that stand for the "normal" or 
standard solar abundance ratio. We notice in Table [T3] 
that the brightening effect is only marginally dependent 
on T e ff and log g. 

We notice that the BCD parameters, D, derived from 
the fitted spectra are consistent with those measured in 
the observed spectra. We can then assert that the ac- 
tual Balmer jump of the studied He-strong stars is well- 
represented by the empirical D parameter. Then, enter- 
ing Table [T2] with the measured D values (Table [7]) in 
columns with high Hc/H ratios, we recover approximately 
the effective temperatures obtained by fitting NLTE mod- 
els (Tabled]). These T off (NLTE) values are similar to those 
reported by other authors (Table [9]). 

b) - Spectral variations 

We observe that HD 37479 (Figure O has near-UV 
flux, Balmer jump and line intensity variations, while the 
Paschen continuum does not seem to undergo appreciable 
changes. The intensity of H lines increases when that of 
He I lines decreases and the near-UV flux is lower. A re- 
markable difference is found in the equivalent widths (~ 
30 %), in the line intensities, and on the He/H line ratio 
as measured in the spectrum obtained on 2004, Sep. 2 
and in another one taken two days before (see Table fTlf . 
In the spectrum of Sep. 2, there is a flux excess in the 
Balmer continuum near the BD, as compared to its level 
in the spectrum of two days latter. This excess is reminis- 
cent of the second component of the BD seen in emission 
in some Be stars (Divan 1978). Whether this flux ex- 
cess in He-strong stars is a matter of an actual emission, 
still remains to be proved. Nevertheless, if we estimate D 
with the criteria explained in § 15.41 a) , we can see that D 
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= 0.094±0.006 dex on 2004, Sep. 2 and D = 0.107±0.004 
dex on 2004, Sep. 4. This implies that the stellar T e g re- 
mains fairly unchanged when there are strong variations 
of He and H lines. Eventhough, the flux excess near the 
BD might lead to too high values of T e ff when derived by 
model fitting. 

HD 37017 displays a similar behaviour as HD 37479 
while HD 64740 does not show appreciable variations for 
the observing dates. 

6 Discussion 

The BCD spectrophotometric classification system has 
shown to be a useful tool to determine fundamental pa- 
rameters for He-weak stars. To this end, we can use the 
calibrations of the BCD (Ai,D) parameters obtained from 
stars near the Sun. 

It is well-known that some He-weak stars display pho- 
tometric and spectroscopic variations and this fact might 
influence the determinations of fundamental parameters. 
However, taking into account that we have observed many 
of the program stars only once, it is interesting to stress the 
good agreement found between the fundamental parame- 
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Table 13: Model monochromatic magnitude brightening 
at A4250 A as a function of He/H, T e g and log g 
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Figure 6: Line variations observed in the spectrum of HD 
37479. The Balmer continuum appears in emission in the 
spectrum taken in Sep. 2, 2004 



ters obtained with the BCD system and with other tech- 
niques corresponding to data taken in different dates and 
probably at different phases. Then we may suggest that, 
in general, the variations due to chemical inhomogenities 
do not introduce significant changes on the BD of He- weak 
stars. However, it would be necessary to study in detail 
those highly variable objects where temperature fluctua- 
tions or variations of the magnetic field might be plaussi- 
ble. Contrary to the He-weak group, He-strong stars show 
large discrepancies among the T e ff values found in the lit- 
erature. There are several arguments that can account 
for these discrepancies: a) methods based on the Kurucz' 
([32)) line blanketed solar composition models can only 
lead to an approximate flux calibration but, if we do not 
have a gravity indicator, it is impossible to choose the tem- 
perature uniquely (Adelman & Pyper p]); b) the use of H/3 
and H7 equivalent widths combined with T c ft calibrations 
and/or theoretical Balmer line fittings are uncertain since 
these lines might be filled by emission or be variable; c) 
Geneva colours are affected by the He-overabundance and, 
furthermore, calibration problems may also be present due 
to the lack of data for hot stars. Zboril et al. ([56]) have 
shown that the [U-B] index of the Geneva photometric sys- 



tem is affected by the enhanced He/H ratio, which affects 
the b-f and b-b He transitions. Due to these difficulties, 
the mentioned methods lead to somewhat uncertain T e g 
determinations . 

In this work we have shown that the BCD calibrations 
lead to overestimated values of T c g for the He-strong stars, 
since they are actually made for objects with He/H w 0.1. 
However, if care is made at identifying the genuine con- 
tinuum flux levels, as indicated in ij5.41 we can still mea- 
sure the BCD parameter D. This parameter not only well- 
represents the actual Balmer jump of He-strong stars but 
also reflects closely the true effective temperature of the 
star and their H/He ratio when it is used in combination 
with results of NLTE models. An alternative way to de- 
termine the effective temperature of He-strong stars could 
be the use of the integrated flux method, where the deter- 
mination of the angular diameter 6 is almost insensitive to 
the type, or quality, of model atmosphere used in the near- 
IR region. The difference between T eff (BCD) and T off (IF) 
translated into BD difference SD may then produce an 
estimate of the He/H ratio in these objects. 

For the He-strong stars, the Balmer jump might be af- 
fected by changes in He/H ratios rather than by changes 
of the photospheric temperature. When increasing the 
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abundance of He, stellar atmospheres of B stars become 
less opaque due to the transparency of He I in the visible 
continuum and the reduction of the relative number of H 
absorbing atoms. Direct calculations of the k + /k~, total 
continuum absorption coefficient ratio at the Balmer dis- 
continuity, show that this effect is slightly stronger for k~ 
than for n + . This leads to smaller values of D for He- 
strong stars than for stars with standard He abundance. 
However, it can still be noted that the He abundance in- 
crease may concern only a fraction / of the apparent stel- 
lar hemisphere. The fit of observed energy distributions 
with models should then use both He/H and / as free 
parameters. The presence of helium patches could be ex- 
plained in terms of the so-called Oblique Rotator Model. 
The enhancement of He/H ratio might also give rise to a 
flux excess in the Balmer continuum originating a second 
component of the BD in emission. 



6.1 Evolutionary status of the studied 
stars 

In order to obtain an insight on the evolutionary status 
of the He-weak and He-strong stars studied in this work, 
we plot in Figure [7] Mboi vs. logT e a. For He-strong 
stars we have taken into account the T c g values derived 
with NLTE models. The isochrones are for normal stars 
of metal content typical for the solar vicinity Z=0.02, and 
were calculated by Bressan et al. ([!]). Some ages i(yr), 
in logarithm scale, are shown in the figure and the curves 
are spaced by steps of 0.1 dex. 

There is a clear distinction made by the effective tem- 
perature, which separates the He-strong from He-weak 
stars. However, most He-weak and He-strong objects are 
in Main Sequence phase. 

Since in Figure [7] there is a well-marked separation 
between He-weak and He-strong stars according to their 
effective temperatures, we are tempted to conclude that 
stars do not evolve from He-weak to He-strong, or vicev- 
ersa. The separation between these stellar classes seems to 
be due mainly to mass-effective temperature related phe- 
nomena. In fact, the known He-variable star HD 125823, 
which is located between the two He-peculiar groups, 
shows alternatively He-w/He-s features. This object could 
provide important clues to disentangle the evolutionary 
status of the He-peculiar stars. 

Our sample of objects is too small to derive correla- 
tions between the star magnetic field and its fundamental 
parameters. However our determinations of fundamen- 
tal parameters and ages are in very good agreement with 
those estimated by Kochukhov and Bagnulo ( [3D] ) • These 
authors find that the surface magnetic flux of CP stars 
increases with stellar age and mass, and correlates with 
the rotation period. 




4.5 4.4 4.3 4.2 4.1 4 



log T afl 

Figure 7: HR-diagram for He-weak stars (dots) according 
to BCD fundamental parameters. We also include the He- 
strong objects (squares) taking into account, in this case, 
the T c ff derived from NLTE models. Isochrones are from 
Bressan et al.([3]). The ages range from 10 6 6 to 10 10 yr 
with steps of A log t = 0.1. 

7 Conclusions 

Several reviews and papers have highlighted the difficulties 
of obtaining the fundamental parameters of CP stars. Of- 
ten, the determination of effective temperatures and sur- 
face gravities of CP stars by means of photometry and 
model fittings are more complex than for normal stars. 

The main purpose of this work is to estimate the fun- 
damental parameters of 28 helium-peculiar stars using 
the BCD spectrophotometric system. The advantage of 
the BCD system is that it provides simultaneously spec- 
tral types, T c ff, logg and bolometric and absolute mag- 
nitudes. All these quantities are obtained from only two 
parameters drawn from low-resolution spectra which are 
not affected by ISM and circumstellar absorptions, and 
the techniques at obtaining them avoid the spectral re- 
gions that are strongly perturbed by chemical anomalies of 
He-peculiar stars. Besides, the effective temperature and 
surface gravity thus obtained are related to deeper atmo- 
spheric layers than those where, on average, the spectral 
lines are formed, since the BCD quantities are obtained 
from the visible continuum spectrum. 

We confirm that the BCD spectrophotometric system 
gives reliable effective temperatures, surface gravities and 
absolute magnitudes for He-weak stars. Not only are T e g 
values in very good agreement with those obtained with 
the UBV, multicolour and Geneva photometry, but also 
with the IRFM and LTE line-blanketed models. Sur- 
face gravities show an excellent agreement with values 
derived from HIPPARCOS parallaxes and stellar evolution- 
ary tracks. The absolute magnitudes also agree with the 
values obtained from the HIPPARCOS parallaxes and the 
mean discrepancy between both methods is on average ± 
0.3 mag. 
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Particular attention is payed in this paper to the effec- 
tive temperature of He-strong stars. The He-strong stars 
are He-variable and this variability not only affects the 
intensity of He and H lines, but also the apparent dis- 
tribution of continuum fluxes. These changes can origi- 
nate difficulties of determining the effective temperatures 
of He-strong stars that were discussed in this work. The 
determination of the BCD parameter D, which quantifies 
the Balmer jump, avoids these inconveniences. Neverthe- 
less, the parameter D is a function of the He/H abundance 
ratio. To specify this dependence, we have calculated non- 
LTE model atmospheres. We have obtained thus a relation 
D = D (He/H) and could determine the effective temper- 
ature of the studied He-strong stars by model fitting. The 
results show that the higher the He/H ratio, the smaller D. 
From these models it is also apparent that increased values 
of He/H are accompanied by brighter fluxes in the Paschen 
continuum. We have also found that the model fitting pro- 
duces much lower effective temperatures than expected 
from the T e s = T c ff(Ai,D) BCD calibrations, which actu- 
ally suit to stars with abundance ratios He/H ks 0.1. We 
have then proposed a method of determining the He/H 
of He-strong stars, which is based on the mentioned D = 
D(Hc/H) relation and on two distinct T g- determinations, 
one coming from the up to now used (Ai,D) calibration, 
and the other, which could be the integrated-flux method. 

We have observed the He-strong star HD 37479, which 
underwent a rather strong variation in the visible spec- 
trum in an interval of 2 days. Although the fitting proce- 
dure of classical stellar atmosphere models with the ob- 
served energy distributions would suggest an apparent 
variation of the effective temperature, a careful determi- 
nation of the D parameter shows that this quantity does 
not vary significantly. We note that in the BCD system D 
can be constant over a large range of Ai parameters, which 
may reflect the fact that the observed spectrophotometric 
variations correspond to changes in the structure of pos- 
sible exophotospheric layers, as in Be stars. In fact, Be 
stars can sometimes exhibit spectrophotometric changes 
implying constancy of the total absolute magnitude My 
(Moujtahid et al. 1998). 

Since the brightening of the Paschen continuum is on 
average not stronger than some AM4250 ~ —0.10 mag, 
the BCD calibration can still be a useful tool to obtain 
a rough estimate of the absolute magnitudes of He-strong 
stars. 

A first insight on the evolutionary status of He- weak and 
He-strong stars was obtained. Both types of He-peculiar 
stars seem to be in the Main Sequence evolutionary phase. 
The He-strong stars are situated roughly in the T c g ~ 
19000 K region of the HR diagram, while the He-weak 
are in the T e g ~ 19000 K zone of the HR diagram for B- 
type stars. More precise Mboi determinations for a large 
number of stars are needed to give a deeper insight on the 
evolutionary status. 

Due to the variability observed in the He-strong group, 



equivalent widths of hydrogen and helium lines are given. 
Data for the He-strong star HD 66765 are reported for the 
first time. Although we have obtained a "normal" Hc/H 
ratio for this star, we should consider that this object was 
reported as a variable He-strong star by Wiegert & Gar- 
rison ([55]). 

Appendix A Line equivalent widths of the 
He-strong stars 

Having in mind the strong variations observed in HD 
37479 and HD 37017, we considered useful to report line 
equivalent width measurements for all the program He- 
strong stars. 

Table 14: Equivalent widths in A of hydrogen and helium 
lines for the program He-strong stars and the Helocentric 
Julian Date 

HD HJD H 7 He 1 
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We measured equivalent widths (EW) of hydrogen and 
helium lines in the spectra of the 8 observed He-strong 
stars and of the variable He-w/He-s star HD 125823 (see 
Table [14")) . We have compared our equivalent width mea- 
surements with those found in the literature. From this 
comparison it comes out that the EWs derived by Wal- 
born ( |54j ) are, in most cases, larger than ours by factors 
ranging from 1.3 to 2, except for the Hei lines observed in 
HD 37479 (Sep. 4, 2004) and HD 264111. The discrepan- 
cies between Zboril's et al. ([5S]) EWs and ours are about 
10 to 15 %. The largest difference is found for HD 37479 
(Sep. 2, 2004), our EWs are larger by a factor of 1.4. 



Acknowledgements: This work was partially sup- 
ported by the Agenda de Promotion Cientffica y Tec- 
nologica (BID 1728 OC/AR PICT 03-12720) and the Pro- 
grama de Incentivos G 11/073 of the National University of 
La Plata. Y. F. acknowledges financial support from the 
Belgian Federal Science Policy (projects IAP P5/36 and 
MO/33/018). This work has made use of CDS data base. 
We would like to thank our referee Dr. Mathys for his 



14 



numerous comments and suggestions which considerably [23] Hubeny, I., & Lanz, T. 1995, ApJ, 439, 875 
helped to improve our manuscript. 

[24] Hubrig, S., North, P., Mathys, G. 2000, ApJ, 539, 352 

References ^ Hun s er ' K -> & Groote, d. 1999, a&a, 351, 554 

[261 Jaschek, C, & Jaschek, M. 1987a, "The Classification 
Adelman, S. J., & Pyper, D. M. 1985, A&AS, 62, 279 ^ Stars „ ; Cambridge University Press , p .l 7 3 

Barbier, D, & Chalongc, D. 1941, Ann. Astrophys., 4, [2?] Jaschek C & Jasche k, M. 1987b, A&A, 171, 380 
30 

[281 Johnson, H.L., & Mitchell, R.I. 1975, Rev. Mexicana 
Blackwell, D. E., & Shallis, M. J. 1977, MNRAS, 180, Astron Astrof 1 299 
177 ' 

[29] Kaufmann, J. P., & Theil, U. 1980, A&As, 41, 271 

Bressan A., Fagotto F., Bertelli G., & Chiosi C. 1993, 
A&AS, 100, 647 [30] Kochukhov, O., & Bagnulo, S. 2006, A&A, 450, 763 

Bychkov, V. D., Bychkova, L. V., & Madej, J., 2003, [31] Kroll R., 1987, A&A, 181, 315 
A&A, 407, 631 [32] KumcZ; R 1979; ApJg; 40; 1 

Cayrel de Strobcl, G., Hauck, B., Francois, P., et al. rool T _ _ „„ „„,,,. .„ „ ... 

1992 A&AS 95 273 [33] Kurucz, R. 1992, CD-ROM No. 19, 20, 21, Cambridge 

Mass.: Smithsonian Astrophysisical Obsrevatoryx 
^Chalongc, D., & Divan, L. 1952, Ann. Astrophys., 15, [M] ^ T ^ A&A ^ ^ m 

™ , ^ o ^. t An a i~> r\ [35] LeBlanc, F., Michaud, G., & Babel, J. 1994, ApJ, 

Chalonge, D., & Divan, L. 1973,A&A, 23, 69 431 388 



Chalongc, D, & Divan, L. 1977, A&A, 55, 117 [36] Leckron6; D ^ FowlCT; j w & Adelman, S. J. 

10] Cidale, L. S., Zorcc J., & Tringaniello L. 2001, A&A, 1974 > A&A ' 32 > 237 



368,160 



[37] Leone, F., & Manfre, M. 1997, A&A, 320, 257 



11] Dimitrijevic, M. S., & Sahal-Brcchot, S. 1984, [3g] Leon6j p ^ CatalanG; F a., & Malaroda, S. 1997, 
JQSRT, 31, 301 A&A; 325; n25 

12] Dimitrijevic, M. S., & Sahal-Brcchot, S. 1990, A&AS, rggn Meg essier, C. 1988a, A&A, 206, 74 
82, 519 

[40] Megessier, C. 1988b, A&AS, 72,551 

13] Divan, L. 1979, in "Spectral Classification of the 
Future", IAU Coloquium 47, eds. Mc Carthy, Philip, [41] Moon, T. T., & Dworetsky, M. M. 1985, MNRAS, 
Coyne, Vatican Observatory, p. 247 217, 305 

14] Divan, L., & Zorec, J. 1982, ESA-SP 177, 101 [42] Moujtahid, A., Zorec, J., Hubert, A. M. et al. 1998, 

A&AS, 129, 289 

15] Flower, P. J. 1977, A&A, 54, 31 

[43] Napiwotzki, R., Schonberner, D., & Wenske, V. 1993, 
16] Freudenstein, S. A., & Cooper, J. 1978, ApJ, 224, A&A, 268, 653 
1079. 

[44] North, P., Jaschek, C, Hauck, B., et al. 1997, Pro- 
17] Gerbaldi, M., Hauck, B., & Morguleff, N. 1974, A&A, ceedings of the ESA Symposium 'Hipparcos - Venice 
30, 105 '97' ; Venice, Italy, ESA SP-402 (July 1997), p. 239-244 

18] Glagolevskij, Y. V. 2002, Bull. Spec. Astrophys. Obs., [45] Osmcr, P. S., & Peterson, D. M 1974, ApJ, 187, 117 
53 33 

[46] Proffitt, C. R., Bragc, T., Lcckronc, D., et al. 1999, 
19] Golay, M. 1972, Vistas in Astronomy, 14, 13 ApJ, 512, 942 

20] Gomez, A. E., Luri, X., Grcnier, S., et al. 1998, A&A, [47] Renson, P., Gerbaldi, M., & Catalano, F. A. 1991, 
336, 953 A&AS, 89, 429 

21] Hauck, B., & North, P. 1993, A&A, 269, 403 [48] Shore, S. N., & Brown, D. N. 1987, A&A, 184, 219 

22] Hauck, B., & North, P. 1982, A&A, 114, 23 [49] Sokolov, N. A. 1998, A&AS, 130, 215 



15 



[50] Sokolov, N. A. 1995, A&AS, 110, 553 

[51] Stepicn, K., & Dominiczak, R. 1989, A&A, 219, 197 

[52] Theodossiou, E., & Danczis, E. 1991, A&AS, 183, 91 

[53] Veto, B., Hempelmann, A., Schoneich, W. & 
Stahlberg, J. 1991, Astron. Nachr., 312, 133 

[54] Walborn, N. 1983, ApJ, 268, 195 

[55] Wicgcrt, P., & Garrison, R. F. 1998, JRASC, 92, 134 

[56] Zboril, M., North, P., Glagolevskij , Y. V., & Bctrix, 
F. 1997, A&A, 324, 949 

[57] Zboril, M. & North, P. 1999, A&A, 345, 244 

[58] Zorec, J., Frcmat, Y., & Cidalc, L. S., 2005, A&A, 
441, 235 

[59] Zorec, J. 1986, PhD. Thesis, Universite de Paris, 
France 

[60] Zorec, J., & Briot, D. 1991, A&A, 245, 150 



16 



